Electronic structure, magnetism and superconductivity of MgCNi 3 
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The electronic structure of the newly discovered supercon- 
ducting perovskite MgCNi3 is calculated using the LMTO and 
KKR methods. The states near the Fermi energy are found 
to be dominated by Ni-d. The Stoner factor is low while the 
electron-phonon coupling constant is estimated to be about 
0.7, which suggests that the material is a conventional type of 
superconductor where Tc is not affected by magnetic interac- 
tions. However, the proximity of the Fermi energy to a large 
peak in the density of states in conjunction with the reported 
non-stoichiometry of the compound, has consequences for the 
stability of the results. 



The perovskite oxides have hosted a wide vari- 
ety of exotic phenomena, including high-temperature 
superconductivity!!, colossal magnetoresistancea and 
ferroelectricitytl. However, the recent report of super- 
conductivity below 8K in the intermetallic compound 
MgCNi3 by He et als is the first in a material which 
has the perovskite structure without any oxygen. Pre- 
sented as a three.-dimensional analog of the borocarbide 
superconductorstj (where there has recently been sugges- 
tion of non-s-wave pairingp), He et al. speculated that 
the favoring of a superconducting ground state over a 
ferromagnetic one in MgCNi3 (with its high Ni content), 
makes it a potential candidate for unconventional super- 
conductivity. In the rare-earth nickel borocarbides, the 
presence of moment-bearing rare-earth atoms and high 
Ni content makes the very existence of superconductivity 
surprising, and some of the phenomena associated with 
its interplay with magnetism have net, been observed Hp. 
any other superconducting materiaC'cl. Huang et al.u, 
however, report no magnetic or structural transitions 
forMgCNi3 in the range [ 2j — 295K. Recent tunneling mea- 
surements by Mao et alS3 have indicated that MgCNi3 is 
a strong-coupling superconductor, and moreover that the 
pairing symmetry could be non-s-wave. Here we present 
calculations of the electronic structure of MgCNi 3 , focus- 
ing on the central issues of magnetism and superconduc- 
tivity. 

The structure of MgCNi3 has been reported to be that 
of a classic perovskitco, comprising a C atom at the body- 
center position, surrounded by a cage of Ni atoms at the 
face-center positions, with the Mg occupying the cube 
corners. The electronic structure of MgCNi3 was calcu- 
lated using both the LMTOEl and, as a cross-check and 



comparison, with the KKPjIJ method. The lattice pa- 
rameter was fixed at the ejtperimental low-temperature 
value of 7.2 a.u. (3.8lA)B. For the LMTO calcula- 
tions, made within the atomic sphere approximation 
(ASA), self-consistency was reached using 286 k-points 
within the irreducible wedge of the simple cubic Brillouin 
zone (BZ). The KKR code used non-overlapping muffin- 
tin potentials, and the mesh of k-points was adaptive 
and determined principally by the required integration 
toleranceaL3. 

The bands, density of states (DOS) and Fermi surface 
(FS) from the LMTO results are shown in Figs. |l|, || and 
|| respectively. In Fig. ||, the dashed line represents the 
contribution of Ni to the total (solid line) DOS. It should 
be noted that near the Fermi energy, Ep, the DOS is 
almost completely due to Ni (see Table |). The KKR 
results were very similar to those from the LMTO, but 
the Ni d-bands were flatter, generating a larger DOS at 
the Fermi level (86.4 compared to 47.6 states (Ry cell) -1 
from the LMTO). Clearly, with the Fermi level lying in 
such close proximity to a large peak in the DOS, there 
will be a sensitivity in the derived quantities. Should 
single crystal samples become available, a determination 
of the FS topology would be a stringent check on the 
location of the Fermi energy on the DOS peak, since the 
low dispersion of the Ni-d bands would make the topology 
very sensitive to small shifts in Ep. 

Henceforth, we will focus on the LMTO results. These 
calculations include s,p,d and f states for all atoms be- 
cause the partial f-DOS is required in the evaluation of 
the electron-phonon coupling. Some bands were found to 
be sensitive to the choice of linearization energies, but a 
separate LMTO calculation using only s-,p- and d-states 
for Ni, and s and p states for Mg and C, gave the same FS 
topology and a very similar total DOS at Ep. Two bands 
cross Ep, making one jungle-gym-like FS sheet around 
the BZ edge, with a spheroid-like sheet around T. The 
second FS plot shows some X-centered shell-like features 
at the BZ faces and small, delicate 'cigars' along T-R. 

The bands (Fig. are quite dispersive near Ep, with 
a Fermi velocity of 2.0 x 10 5 m/s, which in combina- 
tion with the large DOS, would give the material a good 
metallic conductivity. The temperature dependence of 
the resistivity is reported as having the signature of a 
poor metalU As the calculated DOS is strongly vary- 
ing near Ep, it is expected that unusual T-dependencies 
could occur. 

The propensity for a metallic system to adopt a fer- 
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romagnetic ground state can be expressed in terms of 
the Stoner factor, S = 1/(1 — S), this being the ex- 
change enhancement, which diverges at a ferromagnetic 
transitior£3o. The DOS per Ni atom is considerably 
smaller than in fee Ni, and the calculated value of the 
Stoner factor 5=0.43 is far from the ferromagnetic limit. 
This suggests that, in_contrast to common expectations 
of Ni-rich compoundsu, magnetism will not be present 
in this compound and should not interfere much with 
superconductivity. Antiferromagnetic ordering would 
also seem unlikely, since the three dimensional structure 
makes the FS topology quite complicated with no obvious 
nesting features. 

The electron-phonon coupling constant, A, can be ex- 
pressed astj, 



A 



where the sum runs over all atoms, i, with masses, M iy 
and phonon frequencies, u>i, while the numerator, rji, is 
the electronic contribution. Hexe, r\ was calculated in the 
rigid muffin-tin approximationlij. This implies that only 
dipole terms without screening are included. The large 
Ni-d DOS is such that the contribution to r\ is domi- 
nated by Ni d to f scattering. In order to get an idea 
of how large A could be, we estimate the phonon con- 
tribution by using the Debye temperature of fee Ni (450 
kIIZI). In this case, we obtain 0.67 for A, which, with the 
bare DOS, gives an electronic specific heat coefficient of 
4.7 mJ per male Ni K~ 2 . This is smaller than the exper- 
imental values (about 10 mJ per mole Ni K -2 ), which 
leaves room for additional enhancements due to spin- or 
valence-fluctuations. From their specific heat measure- 
ments, He et aln infer A = 0.77 (with errors of +0.17 and 
-0.09). We should like to point out that with the bare 
Fermi energy DOS from the KKR calculations, the agree- 
ment with the electronic specific heat constant would be 
much closer. On the other hand, the LMTO-derived A 
would give a superconducting critical temperature, T c , 
of the order of 10K (using the McMillan formula wi 
/i* = 0.13E3), in line with the experimental observation; 
whereas that derived from the KKR DOS would be much 
larger. It should also be noted that this calculation of A 
is strongly dependent on the phonon frequencies, which 
are not calculated here. 

Since there is a large peak in the DOS close to Ep, 
and if doping (the occupancy of the CUsite is 0.96 rather 
than unity, in the measured materials) were to add 0.5 
holes per unit cell, Ep would reach that peak. In this 
case, the DOS at Ep would be more than twice as large, 
reaching 130 states per cell per Ry, which would mean 
that S would become close to unity. However, this is 
unlikely for the measured material, since the specific heat 
coefficient would become too large, both due to the bare 
DOS and from the additional enhancements. 

The electronic structure of the newly discovered per- 
ovskite superconductor, MgCNi 3 is reported. In sum- 
mary, we find that although the DOS at the Fermi level 



is dominated by the Ni d-states, it is not large enough to 
induce magnetic instabilities. However, the DOS is suffi- 
ciently large to produce strong electron-phonon coupling. 

One of the authors (SBD) would like to thank the 
Royal Society and the DPMC (University of Geneva) for 
financial support. 



1 J. G. Bednorz and K. A. Miiller, Z. Phys. B 64, 189 (1986). 

2 R. M. von Helmolt, J. Wecker, B. Holzapfel, L. Shultz and 
K. Samwer, Phys. Rev. Lett. 71, 2331 (1993) ; S. Jin, T. H. 
Tiefel, M. McCormack, R. A. Fastnacht, R. Ramesh, L. H. 
Chen, Science 264, 413 (1994). 

3 M. E. Lines and A. M. Glass, Principles and Applications 
of Ferroelectrics and Related Materials (Oxford University 
Press, Oxford, 1977). 

4 T. He, Q. Huang, A. P. Ramirez, Y. Wang, K. A. Regan, 
N. Rogado, M. A. Hayward, M. K. Haas, J. S. Slusky, K. 
Inumara, H. W. Zandbergen, N. P. Ong and R. J. Cava, 
Nature (London) 411, 54 (2001). 

5 R. J. Cava et at, Nature (London) 367, 252 (1994) ; R. 
Nagarajan et al, Phys. Rev. Lett. 72 , 274 (1994). 

6 T. Yokoya, T. Kiss, T. Watanabe, S. Shin, M. Nohara, H. 
Takagi and T. Oguchi, Phys. Rev. Lett. 85, 4952 (2000). 

7 P. C. Canfield, P. L. Gammel and D. J. Bishop, Phys. To- 
day 51, 40 (1998). 

8 K. Norgaard, M. R. Eskildsen, N. H. Andersen, J. Jensen, 
P. Hedegard, S. N. Klausen and P. C. Canfield, Phys. Rev. 
Lett. 84, 4982 (2000). 

9 Q. Huang, T. He, K. A. Regan, N. Rogado, M. Hay- 
ward, M. K. Haas, K. Inumaru and R. J. Cava, cond 



mat/0105240 (2001) 



Z. Q. Mao, M. M. Rosario, K. Nelson, K. Wu, I. G. Deac, 
P. Schiffer, Y. Liu , T. He, K. A. Regan and R. G, Cava, 



cond-mat/0105280| (2001) 



11 O. K. Andersen, Phys. Rev. B 12, 3060 (1975) ; T. Jarlborg 
and G. Arbman, J. Phys. F 7, 1635 (1977). 

12 E. Bruno and B. Ginatempo, Phys. Rev. B. 55, 12946 
(1997). 

13 T. Jarlborg and A. J. Freeman, Phys. Rev. B 21, 2332 
(1980). 

14 O. Gunnarson, J. Phys. F 6, 587 (1976). 

15 M. Dacorogna, T. Jarlborg, A. Junod, M. Pelizzone and 
M. Peter, J. Low Temp. Phys. 57, 629 (1984). 

16 W. L. McMillan, Phys. Rev.167, 331 (1968). 

17 C. Kittel, Introduction to Solid State Physics (7 th edition, 
Wiley, New York, 1996). 



2 



0.1 


^-0.1 

& 

W-0.2 
I 

a 

-0.3 
-0.4 
-0.5 




txm r RMR X 

FIG. 1. Electronic band structure of MgCNi3 along the 
high-symmetry directions of the simple cubic BZ. 
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FIG. 2. Density of states (DOS) for MgCNi 3 . The total 
DOS (solid line) and the partial Ni DOS (dashed line) are 
shown. 




FIG. 3. The Fermi surface sheets of MgCNi 3 . 



TABLE I. Site decomposition of the density of states at 
E F [(Ry cell)- 1 ]. 



Site 


s 


P 


d 


f 


Mg 


0.1 


2.0 


0.5 


0.1 


C 


0.2 


3.6 


0.1 


0.2 


Ni 3 


1.0 


2.8 


36.3 


0.3 
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